On the Microscopic Origin of Cholesteric Pitch 
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We present a microscopic analysis of the instability of the nematic phase to chirality when molecular 
chirality is introduced perturbatively. We show that for central force interactions the previously 
neglected short-range biaxial correlations play a crucial role in determining the cholesteric pitch. 
We propose a pseudoscalar strength which quantifies the chirality of a molecule. 
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Chirality in molecules leads to a myriad of macroscopic 
chiral structures, including life itself [Q. A molecule is 
chiral if its mirror image cannot be rotated to replicate 
itself Q|. Equivalently it is chiral if its symmetry group 
does not contain the element S„ - a rotation around an 
Sn axis by 27r/n followed by a mirror through a plane 
perpendicular to that axis. A chiral molecule cannot be 
uniaxial: the only infinite point groups are Coov and Pooh 
1^ which both contain Si. Thus if the molecular orien- 
tations are averaged independently (even if the distri- 
bution is uniaxial about a common molecular axis) the 
interactions will be identical to those of molecules with 
a Coo axis and will therefore not be chiral. Figure |l| 
shows a schematic representation of two interacting chi- 
ral molecules whose degree of chirality can be varied con- 
tinuously as we will discuss below. 




FIG. 1. Schematic representation of two chiral molecules. 
The atoms are represented by both the filled and unfilled cir- 
cles, while the lines serve only as a guide to the eye. Each line 
on each molecule lies in an xz plane parallel to the page. The 
arms and atoms in the plane ai y — L/2 are black and those 
in the plane at y — — L/2 are grey. The angle 7 between the 
projection of the two arms onto the same xz plane determines 
the degree of chirality. As examples, we consider two versions 
of this molecule. In the first, all atoms are identical, while in 
the second, the atoms with a hollow center carry a negative 
charge and those with a filled center carry a positive charge. 
In the nematic phase, the molecules spin freely about the ne- 
matic axis normal to the page so that (sin (j>x) ~ (cos (px) = 0. 
There are, however, orientational correlations between 6 a and 



Molecular chirality induces chiral interactions that pro- 
duce intermolecular torques of a given sign and can 
give rise to equilibrium chiral structures such as the the 



cholesteric phase of liquid crystals. There are two com- 
mon analyses of this effect. The first is purely classical 
iQ-^ , while the second invokes a generalized chiral disper- 
sion force whose origin is quantum mechanical ||^. One 
can argue that there are systems for which the dominant 
interaction is a classical one, involving two-point, central 
forces between atoms or interaction centers on molecules, 
and accordingly in this paper we consider only this classi- 
cal mechanism. In a more detailed paper |^ , we will con- 
sider the quantum interaction and compare its strength 
with the classical one studied here. Our primary focus 
will be the calculation of cholesteric pitch from these in- 
teractions. The usual classical picture of the origin of 
intermolecular twist considers two screw-like molecules 
with excluded volume interactions In order for the 

vanes of the screws to interleave, the molecules must have 
a non-zero angle between their screw axes. A similar pic- 
ture arises via the tangent-tangent interactions of chiral 
molecules |^,7| or via surface-nematic interactions of chi- 
ral dopants |8|. 

The above mechanism produces a preferred nonzero 
rotation angle between long-axes of neighboring chiral 
molecules (such as two screws) only if there is correlation 
between the directions of their short axes, e.g. between 
and (jjg in Figure |^. To our knowledge, no previous 
attempt to calculate the cholesteric pitch has treated this 
properly. If neighboring molecules are spun freely about 
their long axes, they become effectively nonchiral, and 
interactions favoring twist are washed out [|o|. Thus, 
the pitch of a cholesteric depends critically on the degree 
of intermolecular correlation of short axes directions: a 
vanishing correlation leads to an infinite pitch, and com- 
plete correlation, as would be produced by long-range 
biaxial order, leads to the shortest pitch. In a uniax- 
ial phase, mean-field theory does not treat these short- 
axes correlations and cannot predict a finite cholesteric 
pitch from molecular shape. Thus it will either lead to 
a phase with long-range order in the short axes direc- 
tions (a biaxial phase) or to a uniaxial phase without 
chirality Although the results presented here are 

only methodological, they have crucial implications for 
numerical simulations since such calculations are of- 
ten based on excluded- volume hard-core or other central 
classical interactions of the type we consider. In what fol- 
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lows, we will focus on molecules with biaxial symmetry 
and biaxial correlations, though our results can be gen- 
eralized to molecules with n-axial symmetry. Cholesteric 
phases are, of course, biaxial. We are not discussing the 
higher order effect due to the biaxial order induced by 
the cholesteric pitch axis, but rather, the reverse mech- 
anism in which biaxial order is needed to produce chiral 
phases. 

The physics of many chiral liquid crystalline phases can 
be captured via the phenomenological Frank free energy 
density for the unit vector n: 



+ iX3 [n X (V X n)] V ftn- (V X n) , 



(1) 



where Ki , K2 and K3 are the Frank elastic constants for 
splay, twist and bend, respectively, and h is a chiral pa- 
rameter. This free energy is invariant under the inversion 
n ^ — n, consistent with the symmetry of the nematic 
phase. The calculation of h is the focus of this paper. 

The most common manifestation of the preferred pack- 
ing angle of chiral molecules is the cholesteric liquid crys- 
tal phase in which a particular axis of each molecule lies 
along the nematic director n = [cosA;2,sinfcz,0]. The 
z axis is the pitch axis and the pitch is P = 27r/fc, the 
distance over which the nematic director rotates by 27r. 
In this uniform twisted state n • V x n — —k, and the 
Frank free energy density reduces to / = ■^K2k^ — hk 
so that h = —df/dk\k=o- The equilibrium value of k is 
ko = h/K2- Using standard statistical mechanical proce- 
dures, the chiral parameter can be expressed as 
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(2) 



fc=0 



where the brackets denote a thermodynamic average in 
a nematic state in which there is a spatially uniform di- 
rector n, T is the sample volume, and U is the total 
potential energy. We denote the center-of-mass coordi- 
nate of molecule A by and the coordinate relative to 
the center of mass of particle a in molecule A by r^a- In 
general, a should run over all interaction centers (usually 
atoms and nuclei) in molecule A. The potential energy 
U is a function of all the coordinates. To determine h, 
we introduce an infinitesimal twist in a nematic state in 
which the molecules are aligned along a uniform direc- 
tor n. Under such a twist, atomic coordinates within 



molecule A will undergo a rotation Sr], 



where = ke^(e ■ R^) is a rotation angle about an 
arbitrary unit vector e perpendicular to n. The magni- 
tude of Soja increases linearly with the projection of the 
center-of-mass position R^ along e. Using (||) and the in- 
variance of the system with respect to arbitrary rotations 
about the n-axis, we obtain h — — ^T^^ J2ba "^ba, where 
Tba = {'R±-Tba) is the projected torque, Tba is the torque 
exerted on molecule B by molecule A, R = R^ — Ra, and 



Tl± is the projection of R onto the plane perpendicular 
to n. 

The chiral parameter h must be zero in achiral sys- 
tems that have an equilibrium nematic phase. To see 
that our formulation leads to this result, consider eval- 
uating '^^Tba for a system in which all molecules are 
achiral. The thermodynamic average is carried out over 
all configurations of the molecules consistent with the 
assumed nematic order. If, and only if, the molecules are 
achiral, the average may equivalently be carried out over 
all configurations 17, where is obtained from Jl by a re- 
flection through a plane perpendicular to n. But Rj^-x^a 
in ft is the negative of R_l- Tba in fl, so h — ~h = 0. 

Our expression for h is perfectly general: it applies to 
quantum as well as classical systems. Note that hard- 
core interactions can be viewed as the limiting case of 
central forces between atoms of a single kind that mu- 
tually interact via central forces, and thus we begin our 
analysis for such systems. The projected torque is then 



(3) 
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where — r^a and rp — r^p are, respectively, the co- 
ordinates of atoms a and (3 in molecules A and B. To 
facilitate our analysis, we will now consider an expansion 
of Tba in powers of relative atomic distance over center- 
of-mass separation, i.e., in rp,a/R- We expect, however, 
that the conclusions we draw from this analysis are more 
generally valid and apply, in particular, to hard-core in- 
teractions. In such an expansion, only terms that are 
odd in the atomic coordinates and rp are sensitive to 
reflections and thus to chirality. Furthermore, terms that 
are even in and rp are necessarily odd in R and will, 
therefore, not survive the average over the nematic dis- 
tribution function. Thus, if we assume an achiral distri- 
bution characteristic of a nematic phase, we can restrict 
our attention to terms odd in and rp. (We will, how- 
ever, reconsider this point later). Since rp measures the 
position relative to the center of mass, rp = 0, and 
thus the linear and third order terms vanish. The first 
non-vanishing chiral term is the fifth order term 
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(4) 



n (E^^.'^^l''^'';3ar™.rP„r^„9'=a'9™9^'9^y(R) ) , 
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where rpa = rp — ra- This quantity can be reexpressed in 
terms of the second and third rank mass moment tensors 
on molecules X ^ A,B, M^' = J2xex (^{^[-^^l^^''^ 



'jlm 
X 



E 



xex 'x'x'x 



is the usual quadrupole ten- 



sor describing nematic order, which can be decomposed 
into a uniaxial and a biaxial part. If we assume that there 
is perfect alignment of the longest principal axis along the 
nematic direction, then IVP' = Mi{n^n' - i^-'O + Bi\ 
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where i?-'' is the symmetric traceless biaxial tensor with 
no components along n. The third rank tensor S^^^ is 
symmetric. In general, there are correlations between 
the direction of the vector R connecting two molecules 
and the respective orientations of these molecules. The 
important effects we consider are present even if these 
correlations are absent, and we will ignore them. This 
permits us to evaluate the orientational average of prod- 
ucts of W with respect to a distribution that is isotropic 
in the plane perpendicular to n. for example {R^W) — 
Rf^n'n^ + ^Rl{d'^ -n'n^). Setting ^(R) = g{R^/2), ex- 
panding the derivatives of V , and performing the above 
average, we obtain 



TL = b^,kQ'' {{B^Jsr + B^jSl^') X(R)) , (5) 
where Q^^ = ji^-nP — and 

if (R) = i?i{.g(^) + + [.gW + Rlg^'^] }, (6) 

where g'^"'\x) — d^g{x)/dx'^. We see then, that only the 
traceless part of S'^^'' contributes to and so we may 
take it to be traceless. 

As we have already discussed, T^a is nonzero only if 
the molecules are chiral. How is this fact manifested 
in (H)? Since both Q^p and S-'' can be nonzero for 
achiral molecules, it would seem that the tensor S''^'' 
is a measure of chirality. This is not true, however, 
because S^^^ also has components that can be nonzero 
for achiral molecules. Though there are many possi- 
ble definitions of a molecular chiral strength, when a 
molecule has a unique long axis we propose the pscu- 
doscalar ijj — S^^"^eijkQ^^ B^™" as a measure of the chiral 
strength of a molecule. This definition is useful because, 
as we shall see, h is proportional to 'tp. If the molecule 
is not biaxial {i.e., B^™ = 0), V will vanish. In addi- 
tion, since '0 is a rotational invariant odd in r, it will 
also vanish if the molecule is achiral. It is possible that 
even for a chiral molecule, '0 vanishes. If this were the 
case, however, there would still be nonvanishing contri- 
butions to at higher order in the expansion in powers 
of r / R. Indeed, a complete description of chiral interac- 
tions requires the knowledge of all the chiral moments of 
the molecules. In the basis of the principal axes of the 
molecule with n along z,ip = S''t"'{B'='' - Ryy). Only the 
components of S''^'^' with three different indices in this ba- 
sis contribute to 0. We can, therefore, replace S''^'' in (|^) 

with the tensor S whose only nonvanishing component 
in the principal axis basis is S^^^. 

The projected torque Tj^ is an average over fluctu- 
ations in the aligned nematic phase. It will be zero 
unless biaxial directions on pairs of molecules are cor- 
related. If we were to spin the molecules indepen- 
dently, then (B^™) would be zero, and both terms in 
would vanish. However, when there are biaxial corre- 
lations, this term will not vanish. To see this explicitly. 



we can use the identity s'^'' = —ipA^prsB'^^Q^'' /2B'^ + 
(5 symmetric permutations), where B^ = Tr(i3^). In- 
serting this into (||), we come to (assuming identical 
molecules) : 

= ^ {KmB^^B'^) EE ^ (if(R)rb(R)) (7) 

where rb(R) oc (cos [2 (0^ — 4>a)]) is the biaxial correla- 
tion function for two molecules separated by a distance R. 
We would expect in a uniaxial phase that T}){R) cx e~^^^ 
where ^ is the biaxial correlation length. Naively, one 
would expect ^ to be of the order of the molecular spac- 
ing. Thus we conclude that at the very least chirality 
requires biaxial correlations among the nematogens. We 
tabulate for a number of chiral molecules in Table |. 

TABLE I. Value of for molecules made of atoms located 
at the coordinates given in the first column. The first molecule 
is shown in Fig. 1 and the second is a helix of uniform density. 



Atomic Coordinates 



L>2w: {{±w,0,-L/2), 

{itw cos 7, ±w sin 7, L/2)} 

L->r,neZ: | s G [-1/2, 1/2], 
(j- cos(27rns), r sin(27rns), Ls^ | , 



2w''Lsin(27) 
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(7rn) 



Equation (^ gives the dominant contribution to Tg^ to 
linear order in tp in the nematic phase. There is an addi- 
tional contribution linear in -0 arising from the chiral part 
of the equilibrium probability distribution {e~^ 1^^^ j Z 
where Z is the partition function) and those terms with 
even powers of r arising in r^OkV that are averaged in 
(^. In the isotropic phase, our analysis can be extended 
to show that this term cancels the contribution to Tg^ 
from (^) and (^ to produce ft, = as required. In the or- 
dered phase, this other term is higher order in correlation 
functions and is subdominant to our result 

We have shown that the projected torque Tg^ and 
hence K2kQ will be proportional to the molecular chiral 
strength ip. We note that there are a number of chiral 
liquid crystals, such as solutions of the viruses FD and 
TMV as well as of DNA, that show very small, if any, 
macroscopic chirality ||l^,0. We believe the ideas pre- 
sented here explain these observations, although a com- 
plete understanding will require a thorough investigation 
of the quantum dispersion force. Helical molecules have 
very small biaxiality and hence small values of (see Ta- 
ble |) . The chiral contribution to ^p comes from S'^^^ and 
is inversely proportional to the number of turns per unit 
length, which is consistent with one's geometric intuition. 
In addition, since the molecular chiral strength depends 
on the degree of molecular biaxiality, we see that for fixed 
turn density (or equivalently, fixed S^p''), ip falls off as 
the total length squared. In addition, these molecules 
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can easily rotate independently and are far apart (tens 
of Angstroms). Hence we believe the biaxial correlation 
length ^ will be small compared to the intermolecular 
spacing. Our classical analysis should be valid especially 
for FD and TMV since these molecules are thought to 
interact sterically. Alternatively, short-molecule, ther- 
motropic liquid crystals show very strong chirality, with 
pitches on the order of 5000A. These molecules are gener- 
ally quite flat and thus quite biaxial. Typical molecular 
densities would not allow the molecules to rotate inde- 
pendently of each other, and thus we expect the biaxial 
correlations at the molecular separation to be reason- 
ably large. The combined effect of a large biaxial com- 
ponent to tp and of a large ^ should lead to relatively 
short pitches. In both cases we note that naively one 
would expect on dimensional grounds that kg would be 
on the order of 7r/a where a is the intermolecular spacing, 
which is certainly not a typical inverse pitch. Pitches are 
typically on the scale of microns, not Angstroms. Our 
expression (0) for the leading term of Tg^ is consistent 
with all of the above observations. We also note that in 
all but the most dilute solutions we do not expect any 
universal dependence of pitch on concentration or tem- 
perature: the details of the interactions and correlations 
should be different from system to system. 

We briefly mention a number of generalizations to be 
discussed later ||^. We have considered here only the 
interactions between a pair of perfectly aligned, identical 
molecules. In the nematic state, the molecules are not 
perfectly aligned and the Maier-Saupe order parameter 
S is less than 1. We can incorporate these fluctuations 
into the calculation of T^a. Indeed, we find, as discussed 
above, that when there is no nematic order there is no 
net torque. Since (|^) involves the product of S'^'p^ on 
one molecule and B^' on the other, our results easily 
generalize to chiral molecules interacting with achiral, 
biaxial molecules. More generally, we find that including 
correlations between the intermolecular direction and the 
molecular orientation leads to chiral interactions between 
chiral molecules and uniaxial molecules. 

Additionally, we note that atomic identity may be 
relaxed. In this case, there is an interaction Vc be- 
tween pairs of atoms leading to a potential energy U = 
iZ]9q9/3^c(R- + i"b/3 - r^c), where Qa and qp are the 
"charges" of atoms a and (3. A chiral molecule such 
as that shown in Figure ^ can carry a dipole moment 
Pa = X]q la'^^Aa perpendicular to its long axis. We 
then find a third order contribution to the projected 
torque of the form eijfep^Ci™i?^9'^9'9™V^, where C-'^ — 

E/3 9/3 ('^^^/T ~ i^^^"*™) charge quadrupole moment 

tensor. The molecular chiral strength analogous to 
for this system is V'c = eijkP^C^"^Q^^ . With this defi- 
nition, Tba is proportional to tpc^p{R)) where rp(R) is 
the dipole-moment pair correlation function, which, like 
rf,(R) measures angular correlations in the plane perpen- 



dicular to n. Adding charges q and —q to the molecules 
shown in Figure [l](and using the atomic positions spec- 
ified in Table |), we have V'c = 4q^w^Lsin7. We believe 
that chiral interactions of this type play a role in those 
liquid crystals which can form ferroelectric phases (i.e., 
Sm-C*, TGB^, etc.). 

We close with some observations concerning the re- 
lation of our work to previous treatments of chiral in- 
teractions. An intermolecular potential of the form 
= M'JejkiR''Ml%{R) + {A ^ B) remains chiral 
and nonzero upon spinning about the local nematic di- 
rector. It leads automatically to a free energy of the 
form of (|l|) with h proportional to Vp. Thus, this po- 
tential or ones similar to it are often used as a starting 
point for the description of chiral liquid crystals. Our 
analysis shows that this potential cannot be obtained 
from classical central forces between atoms on molecules, 
though it can arise through quantum dispersion forces 
The potential corresponding to ^ has the form 
VsA = Y.BA[sTM''^d'did''d^d"'V + [A ^ B)]. S'^^ 
is a symmetric tensor - it cannot be expressed in terms 
of tijk and M^'- to produce V^\. 
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